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Abstract 

Background: The association between polymorphisms rs6265 and rs2030324 in brain-derived neurotrophic factor (BDNF) 
and Alzheimer's disease (AD) has been widely reported, but the results remain controversial. 

Methods: A comprehensive search of Pubmed, Web of Science, China National Knowledge Infrastructure (CNKI), Wanfang 
Med Online and China Biology Medical literature database (CBM) was performed. Pooled odds ratios (ORs) with 95% 
confidence intervals (CIs) were calculated using fixed or random-effects models. We excluded the studies with OR>3.0 or 
OR<0.3 for sensitive analysis. Subgroup analysis by ethnicity, form of AD and gender was carried out. Meta-regression was 
conducted to explore the potential sources of between-study heterogeneity. 

Resu/ts: 29 articles with 7548 cases and 7334 controls concerning rs6265 and 22 articles with 5796 cases and 5706 controls 
concerning rs2030324 were included in this meta-analysis. The combined evidence suggested rs6265 contributing 
significantly to the increased risl< of AD in females (codominant: fixed-effects model (FEM): OR = 1.13, 95% CI = 1.04-1.23; 
dominant: FEM: OR = 1 .1 7, 95% CI = 1 .05-1 .31 ), especially for Caucasian females (codominant: FEM: OR = 1 .1 8, 95% CI = 1 .03- 
1.34; dominant: FEM: OR = 1.18, 95% CI = 1.01-1.37) and female late-onset Alzheimer's disease (LOAD) patients (codominant: 
FEM: OR =1.22, 95% CI = 1 .05-1 .41; dominant: FEM: OR =1.23, 95% CI = 1 .03-1 .46). No evidence indicated an association 
between rs2030324 with AD in codominant (random-effects model (REM): OR = 1.06, 95% CI = 0.89-1. 26) and dominant 
(REM: OR = 1.05, 95% CI = 0.86-1 .27) models. 

Conclusion: J\r\\s meta-analysis suggested A allele of rs6265 might increase the risk of AD in Caucasian females and female 
LOAD patients. In addition, no evidence indicated an association between rs2030324 with AD. Further studies are needed to 
confirm these results. 
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Introduction 

Alzheimer's disease (AD) is an age-associated neurodegenerative 
disorder characterized by progressive decline in cognitive function, 
which typically begins with deterioration in memory [1]. The 
number of people with dementia worldwide in 2010 is estimated at 

35.6 million and is projected to nearly double every 20 years to 

65.7 miUion in 2030 and 115.4 miUion in 2050. AD is the most 
common form of dementia and possibly contributes to 60-70% of 
cases [2]. In the US alone, AD is related with an estimated health- 
care cost of US$172 billion per year [3]. The overwhelming 
number of AD patients, combined with the staggering economic 
burden, makes AD a public health problem. 

The key pathological changes that observed in AD brain tissue 
are the accumulation of neuritic extracellular amyloid plaques and 
intracellular neurofibrillary tangles [4]. However, the neuropath- 
ological etiology of AD remains unclear, but are probably related 



with the combined interaction between gene variants and 
environmental factors [5]. There have been many genetic 
polymorphisms reported to be associated with AD, such as 
amyloid precursor protein (APP), presenHin- 1 (PffiA^), presenUin- 
2{PSEN2) [6], apolipoprotein E [APOE) [7,8] and sortilin-related 
receptor 1 (SORLl) [9] . However, as a complex disorder, the genes 
mentioned-above can not explain the overall genetic susceptibility 
and it is supposed that some other genes may participate in the 
development of AD. 

Brain-derived neurotrophic factor (BDNF), as a member of the 
neurotrophic family, plays an important role in the growth, 
development, differentiation and regeneration of various types of 
neurons in the central nervous system [10]. Autopsy studies found 
reduced mRNA expression of BDNF in the hippocampus of 
patients with AD [1 1], which imphcates the possible participation 
of BDNF in the pathogenesis of AD. There have been many 
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polymorphisms studied in BDMF gene, such as rs 11 030 104, 
rsl6917204, rs7103411, rs6265 and rs2030324. However, only 
the last two polymorphisms have been widely studied, with no 
linkage disequilibrium (LD) between them. What's more, these 
results are inconsistent and individual studies have relatively small 
power to confirm this association. For example, the G allele of 
rs6265 confers risk elfect for AD in subjects of Japanese 
(OR= 1.23, 95% CI= 1.02-1.47) [12], but no significant associ- 
ation was found in Italians [13]. Therefore, we performed a meta- 
analysis to identify the association of the two polymorphisms in 
BDNF and AD susceptibility. 

Materials and Methods 

Search Strategy 

A literature search was performed for available articles that 
were published in English or Chinese (up to November 2013) from 
the following databases: (1) Pubmed; (2) Web of Science; (3) China 
National Knowledge Infrastructure (CNKI); (4) Wanfang Med 
Online; (5) China Biology Medical literature database (CBM). The 
search used the following keywords: "Alzheimer's disease" or 
"AD" and "brain-derived neurotrophic factor" or "BDNF' and 
"polymorphism" or "mutation" or "variant". We also reviewed 
the references of included articles to identify additional articles not 
captured by our database searches. 

Inclusion Criteria 

Two investigators reviewed all relevant studies independently to 
determine whether an individual study was eligible for inclusion. If 
the two investigators disagreed about the eligibility of a study, a 
senior researcher was invited to the discussion. The inclusion 
criteria were as follows: (1) case-control or cohort study published 
as an original study to evaluate the association between rs6265 and 
rs2030324 polymorphisms in BDNF gent and AD susceptibility; (2) 
AD were diagnosed according to NINCDS— ADRDA criteria, 
DSM-IV criteria or CERAD criteria; (3) genotype frequencies 
were reported in the articles or could be calculated or ORs and 
95% CIs can be obtained; (4) the genotype frequencies of controls 
are consistent with Hardy- Weinberg equilibrium (HWE); (5) if one 
study from the same population had been published more than 
once, we choose the most recent or complete one; (6) English or 
Chinese language articles were included. 

Data Extraction 

Two investigators extracted the data independently and reached 
a consensus on all items. Information extracted from each study 
was as follows: first author, publication year, country, ethnicity of 
studied population, diagnostic criteria for cases, sample size, 
genotype distributions, mean age and the percentage of male. 

Statistical Analysis 

The chi-square [f^] analysis was used to test deviation from 
Hardy-Weinberg equilibrium (HWE) for the rs6265 and 
rs2030324 genotype distribution of BDNF gene in control groups, 
and P<0.05 was considered as departure from HWE. Pooled 
measure was used as the inverse variance weighted mean of the 
logarithm of odds ratio (OR) with 95% confidence intervals (CI) to 
evaluate the strength of the association of rs6265 and rs2030324 in 
BDNF gene with risk of AD. We conducted analysis for each 
polymorphism considering dominant (AA-I-GA vs. GG for rs6265, 
TT-I-CT vs. CC for rs2030324), recessive (AA vs. GA-I-GG for 
rs6265, TT vs. CT+CC for rs2030324) and codominant (A vs. G 
for rs6265, T vs. C for rs2030324) models, respectively. I^ of 
Higgins and Thompson was used to describe heterogeneity among 



studies [14]. The random- effects model (REM) was adopted if 
significant heterogeneity (I^>50%) was found; otherwise, the 
fixed-effects model (FEM) was adopted, and Mantel-Haenszel was 
used to assess the fixed effects. For rs6265 and rs2030324, meta- 
regression with restricted maximum likelihood estimation was 
performed to explore the potentially important covariates that 
might exert substantial impacts on between-study heterogeneity. 
Specific genetic variants causally associated with common diseases 
would have small effects (risk ratios mostly <2.0) [15,16]. 
Therefore, for sensitive analysis, we excluded the studies with 
OR>3.0 or OR<0.3 for both of the two polymorphisms (rs6265 
and rs2030324) to control the impact of outher values resulting 
from low cell counts within each single study on the pooled effect. 
Moreover, the 'leave one out' sensitive analysis was performed 
using I^>50% as the criteria to evaluate the key studies with 
substantial impact on between-study heterogeneity [17]. When 
heterogeneity was observed, subgroup analysis was also carried 
out. Pubhcation bias was evaluated by Harbord's test. An 
influence analysis was performed to describe how robust the 
pooled estimator is to removal of individual studies. If the point 
estimate of its omitted analysis lies outside the 95% CI of the 
combined analysis, the individual study is suspected of excessive 
influence [18]. All the statistical analyses were conducted using the 
STATA version 10 (Stata Corporation, College Station, TX, 
USA). Two-tailed PsO.05 was considered as statistically signifi- 
cant. 

Results 

Literature Search and Study Characteristics 

The search strategy identified 309 articles from English 
databases (177 articles from Web of Science, 132 articles from 
Pubmed), 176 articles from Chinese databases (161 articles from 
China National Knowledge Infrastructure, 9 articles from 
Wanfang Med Online, 6 articles from China Biology' Medical 
literature database). 64 articles were reviewed in fuU-text. 
Furthermor(;, w(; c-xcludcd 30 artick-s according to the inclusion 
criteria, and obtained 1 additional article through references 
review. Finally, 35 articles were included in this meta-analysis. 
Figure 1 showed the flow diagram of literature search. 

For rs6265, data from 29 published articles [12,13,19-45] with 
32 studies were used including 7548 cases and 7334 controls. For 
rs2030324, 22 articles[12,13,22,25,26,29,31-36,38-40,43,46-51] 
with 23 studies including 5796 cases and 5706 controls were used. 
The ethnicity of the studies included Caucasian, Asian, African 
and others. AH articles were case-control studies. The detailed 
characteristics of the two polymorphisms are showed in Table 1 , 
Table 2, Table 3 and Table SI in File SI. 

Influence Analysis and Publication Bias 

For both of the two polymorphisms, no individual study has 
excessive influence on the pooled effect in any of dominant, 
recessive and codominant models (data not shown). Harbord's test 
showed no publication bias for both two polymorphisms. 

Quantitative Synthesis 

Overall analysis for rs6265 and rs2030324. For overall 
analysis, no significant association was found between A allele and 
risk of AD in codominant (REM: OR= 1.03, 95% CI = 0.95- 
1.12), dominant (REM: OR= 1.06, 95% CI = 0.95-1.18) and 
recessive model (FEM: OR = 1.00, 95% CI = 0.89-1.12) for 
rs6265 polymorphism. No evidence indicated an association 
between rs2030324 with AD in codominant (REM: OR= 1.06, 
95% CI = 0.89-1.26), dominant (REM: OR = 1.05, 95% 
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Figure 1. Flow diagram of literature search. 
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CI = 0.86-1.27) and recessive (FEM: OR= 1.09, 95% CI = 0.8,5- 
1.39) models. The results for overall analysis are showed in 
Table 4. Figure 2 presented the forest plot of ORs in codominant 
model (A vs. G) in overall analysis for rs6265. 

Evidence for heterogeneity (I >50%) was found in codominant 
and dominant models considering the association of rs6265 and 
rs2030324 polymorphisms with AD. Univariate meta-regression 
with the covariates of publication year, ethnicity, diagnostic 
criteria, form of AD, age (ratio of mean age in case group to that in 
control group) and gender (ratio of male percent in case group to 
that in control group) for the above-mentioned polymorphisms 
showed that no covariates have significant effect on between-study 
heterogeneity. 



Subgroup Analysis for Both Two Polymorphisms 

In the stratified analysis by ethnicity, no evidence indicated the 
association between rs6265 and AD susceptibility both for 
Caucasians and Asians in codominant (Caucasian: REM: 
OR= 1.03, 95% CI = 0.90-1.17; Asian: FEM: OR= 1.03, 95% 
CI = 0.95-1.11), dominant (Caucasian: REM: OR = 1.04, 95% 
CI = 0.88-1.22; Asian: FEM: OR= 1.05, 95% CI = 0.94-1.18) 
and recessive (Caucasian: FEM: OR= 1.01, 95% CI = 0.81-1.25; 
Asian: FEM: OR= 1.01, 95% CI = 0.88-1.15) models. After 
stratified by gender, a statistical significant association between 
rs6265 and risk of AD was observed for females in codominant 
(FEM: OR=1.13, 95% CI = 1.04-1.23) and dominant model 
(FEM: OR=1.17, 95% CI= 1.05-1.31), but not in recessive 
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model (FEM: OR=1.13, 95% CI = 0.95-1.35). However, no 
association was detected for males in codominant model (FEM: 
OR=0.98, 95% CI = 0.88-1.08), dominant model (OR=1.01, 
95% CI = 0.88-1. 16) and recessive model (FEM: OR = 0.91, 95% 
CI = 0.74-1.12). When stratified by ethnicity in females, for 
Caucasian females, the A allele was found contributing signifi- 
candy to the increased risk of AD in codominant model (FEM: 
OR =1.18, 95% CI = 1.03-1.34) and dominant model (FEM: 
OR= 1.18, 95% CI = 1.01-1.37), not in recessive model (FEM: 
OR= 1.40, 95% CI = 0.94-2.10). With regard to Asian females, 
no association was detected in codominant model (FEM: 
OR = 1.09, 95% CI = 0.98-1.22), dominant model (FEM: 
OR=1.15, 95% CI = 0.98-1.36) and recessive model (FEM: 
OR= 1.09, 95% CI = 0.90-1.33). When stratified by gender in 
LOAD patients, the A allele was observed significantly associated 
with AD in female LOAD patients in codominant model (FEM: 
OR= 1.22, 95% CI= 1.05-1.41) and dominant model (FEM: 
OR= 1.23, 95% CI = 1.03-1.46), but not in recessive model 
(FEM: OR= 1.47, 95% CI = 0.88-2.44). However, no association 
was found in male LOAD patients in any of the above-mentioned 
models. Figure 3 presented the forest plot of ORs in codominant 
model (A vs. G) in female group for rs6265. 

For rs2030324, in the subgroup analysis, for Caucasian, no 
association was found in codominant model (REM: OR = 1 .00, 
95% CI = 0.82-1.21), dominant model (REM: OR = 0.97, 95% 
CI = 0.78-1.22) and recessive model (FEM: OR=1.10, 95% 
CI = 0.86-1.40). With regard to Asians, there was no association 
between T allele and AD susceptibility in codominant model 
(REM: OR= 1.22, 95% CI = 0.80-1.85), dominant model (REM: 
OR = 1.24, 95% CI = 0.82-1.89) and recessive model (FEM: 
OR= 0.81, 95% CI = 0.18-3.71). When stratified by gender, there 
was still no association between rs2030324 and AD both for 
females and males in codominant model (Female: FEM 
OR=0.97, 95% CI = 0.75-1.26; Male: FEM OR = 1.02, 95% 
y CI = 0.74-1. 40), dominant model (Female: FEM: OR = 1.02,95% 

|3 CI = 0.74-1.40; Male: OR = 0.87, 95% CI = 0.56-1.34). The 

^ pooled results are summarized in Table 5 and Table S2 and 

t Table S3 in File SI. 

"a) 

i Sensitivity Analysis 

.| After excluding articles with OR>3.0 or OR<0.3, low 

J heterogeneity (I^<50%) was found in the codominant, dominant 
and recessive models, and the results for both two polymorphisms 

u are consistent with the ones without sensitive analysis. Further- 

> more, the "leave one out" analysis was carried out when 

t significant heterogeneity (I^>50%) was observed. After excluded 

P key articles with substantial impact on between-study heterogene- 

-§ ity, the data showed rs6265 confers a risk elfect for Caucasians in 

i codominant model (FEM: OR= 1.08, 95% CI = 1.00-1.17, 

I P= 0.042). With respect to rs2030324, the result indicated T 

I allele was associated with EOAD patients in codominant (FEM: 

Q OR=2.02, 95% CI= 1.18-3.47) and dominant model (FEM: 

>^'S OR= 2.01, 95% CI= 1.15-3.52). 

% I Discussion 

o o 
£ °. 

^ ^ To our knowledge, many case-control studies have been carried 

c q- out to investigate the role of BDMF gene in the development of 

o 1 AD. However, these results remained controversial. 29 articles 

Q concerning rs6265 and 22 articles concerning rs2030324 were 

^ ^ included in our meta-analysis. With respect to rs2030324, there 

S 5 was no evidence for an association with AD. What's more, the 

§ ;5 combined evidence suggested that rs6265 was not associated with 
AD for overall analysis. After stratified by ethnicity in females, our 
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Author 



Year 



OR (95% CI) 



Weight {%) 



Caucasian 
Boiocchi et al. [13] 
Borronietal. [21] 
Feher et al. [24] 
Cozza et al. [26] 
Huang etal. [29] 
Zhang et al. [31] 
Saarela et al. [34] 
Desai el al. [36] 
Vepsalainen et al. [38] 
Bodner et al. [40] 
Li et al.[41] 
Li et al. [41] 
Li et al. [41] 
Nacmias et al. [42] 
Bagnoli et al. [43] 
Combarros et al.[44] 
Ventriglia et al. [45] 
D+L Subtotal (l-squared ■ 
l-V Subtotal 

Asian 

Sonali etal. [19] 
Ou et al. [20] 
Fukumoto et al. [22] 
Qi et al. [23] 
Qian et al. [25] 
Yu et al [27] 
He et al. [28] 
Tsai etal. [32] 
Akatsu et al. [33] 
Bian et al. [37] 
Matsushita etal. [12] 
Nishimura et al. [39] 
D+L Subtotal (l-squared ■ 
l-V Subtotal 



2013 
2012 

2009 » 

2008 

2007 

2005 

2006 

2005 

2005 

2005 

2005 

2005 

2005 

2004 

2004 

2004 

2002 

■■ 64.6%, p = 0.000) 



2013 
2012 
2010 
2009 
2008 
2008 
2007 
2006 
2006 
2005 
2005 
2005 
= 29.8%, 



4— 



H4- 



p = 0.154) 



D+L Overall (l-squared = 54.0%, p = 0.000) 
l-V Overall 

NOTE: Weights are from random effects analysis 

.281 



0.97 
1.06 
0.39 
1.10 
1.28 
1.17 
1.90 
1.08 
0.99 
1.04 
1.08 
1.19 
1.15 
1.06 
1.26 
1.11 
0.66 
1.03 
1.03 



1.19 
0.97 
1.17 
1.25 
0.96 
1.08 
0.98 
1.39 
1.00 
1.11 
0.82 
0.89 
1.04 
1.03 

1.03 
1.03 



0.73, 
0.77, 
0.28, 
0.73, 
0.83, 
0.87, 
1.08, 
0.90, 
0.75, 
0.75, 
0.83, 
0.88, 
0.88, 
0.65, 
0.82, 
0.80, 
0.44, 
0.90, 
0.95, 



0.70, 
0.57, 
0.99, 
0.81, 
0.66, 
0.73, 
0.82, 
1.04, 
0.68, 
0.85, 
0.68, 
0.68, 
0.94, 
0.95, 



1.28) 
1.46) 
0.55) 
1.64) 
1.97) 
1.58) 
3.32) 
1.29) 
1.32) 
1.46) 
1.40) 
1.62) 
1.50) 
1.71) 
1.93) 
1.54) 
1.00) 
1.17) 
1.11) 



2.02) 
1.65) 
1.38) 
1.94) 
1.41) 
1.61) 
1.15) 
1.86) 
1.49) 
1.45) 
0.98) 
1.18) 
1.14) 
1.11) 



0.95, 1.13) 
0.97, 1.08) 



3.95 
3.48 
3.32 
2.71 
2.46 
3.76 
1.72 
5.35 
3.98 
3.36 
4.22 
3.70 
4.15 
2.13 
2.51 
3.44 
2.64 
56.88 



1.85 
1.87 
5.57 
2.45 
2.88 
2.78 
5.52 
3.82 
2.79 
4.17 
5.34 
4.05 
43.12 



100.00 



3.56 



Figure 2. Forest plots of relationship between BDNF gene rs6265 polymorphism and AD risk in codominant model (A vs. G) for 
overall analysis. White diamond donates the pooled OR. Black squares indicate the OR in each study, with square sizes inversely proportional to the 
standard error of the OR. Horizontal lines represent 95% CIs. 
doi:1 0.1 371 /journal.pone.0094961 .g002 



data indicated rs6265 lead to the increased risk of AD in 
Caucasian females, but not for Asians. When stratified by gender 
in LOAD patients, the A allele was found contributing significantly 
to the increased risk of AD in female LOAD patients. The results 
of our study were not consistent with the previous meta-analysis 
conducted by Fukumoto et al. [22] in 2010 including 16 studies 
revealing a gender-related association between rs6265 polymor- 
phism and AD susceptibUity. The reason might be that the 
quantitative assessments in our study were based on a larger 
sample size and we have performed a detailed subgroup analysis. 

The Met66-BDNF protein has been shown to be associated 
with reduced transport of BDNF from the Golgi region to 
appropriate secretory granules in neurons, compared with the 



Val66-BDNF protein [52]. Moreover, the A allele of rs6265 was 
related with poorer episodic memory, abnormal hippocampal 
activation, and lower hippocampal n-acetyl aspartate (NAA) in 
human subjects [53]. Epidemiological studies showed higher 
incidence and prevalence of AD in women than in men [3,5]. 
Molecular mechanisms underlying this correlation had been 
considered as colocalization estrogen receptors with BDNF- 
synthesizing neurons in the forebrain [54] and induction BDNF 
expression by estrogen through the estrogen response element 
[55]. In addition, though the two forms of AD (EOAD and 
LOAD) have different patterns of genetic epidemiology [4], the 
results of our study revealed no difference between them. 



PLOS ONE I www.plosone.org 



8 



April 2014 I Volume 9 | Issue 4 | e94961 



BDNF Gene and Alzheimer's Disease 



Author 
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Figure 3. Forest plots of relationship between BDNF gene rs6265 polymorphism and AD risk in codominant model (A vs. G) for 
female group. White diamond donates the pooled OR. Black squares indicate the OR in each study, with square sizes inversely proportional to the 
standard error of the OR. Horizontal lines represent 95% CIs. 
doi:10.1371/journal.pone.0094961.g003 



Between-study heterogeneity is common in meta-analysis for 
genetic association studies [56] and it is essential to explore the 
potential sources of between-study heterogeneity [57]. This meta- 
analysis also showed significant between-study heterogeneity in 
dominant model and codominant model for both two polymor- 
phisms. An indeterminate number of characteristics that varied 
among studies could be the sources of between-study heterogene- 
ity, such as publication year, ethnicity, diagnostic criteria, form of 
AD, age, and gender etc. Therefore, in order to explore the 
potential sources of between-study heterogeneity for rs6265 and 
rs2030324, meta-regression was adopted. However, the afore- 
mentioned covariates were not important contributors to this 
disease-effect heterogeneity. Considering that our meta-analysis 
showed significant heterogeneity, subgroup analyses by ethnicity 
(Caucasian and Asian), form of AD (EGAD and LOAD) and 
gender (female and male) etc. were performed to explore the 
sources of heterogeneity. However, between-study heterogeneity 
still existed in subgroups, suggesting the presence of other 
unknown confounders. AD is a complex multi-factorial disease 
and is related with the combined effects between gene variants and 



environmental factors. Therefore, other genetic and environment 
variables, as well as their possible interaction, may be potential 
contributors to this disease-effect unconformity. We further 
conducted a sensitivity analysis excluding articles with OR>3 or 
OR<0.3. After sensitivity analysis, low heterogeneity (I^<50%) 
was found in the codominant, dominant and recessive models, and 
the results are consistent with the one before heterogeneity analysis 
for rs6265 and rs2030324, strongly identified the stability of our 
results. Moreover, no publication bias was found in any of the 
above-mentioned models for both two polymorphisms. 

The major strength of our meta-analysis is that the results were 
based on the large number of participators, allowing a much greater 
possibility of reaching definitive conclusions. Additionally, we 
conducted subgroup analyses to explore the potential sources of 
heterogeneity, and sensitivity analysis was carried out to ensure the 
stability of our results. However, our meta-analysis also had some 
limitations. Firstiy, lack of the original data of included articles made 
it impracticable to excluded potential confounders completely, 
especially the confounding of age. Secondly, different diagnostic 
criteria may have possible influence on the diagnosis of AD. 
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In conclusion, this meta-analysis suggested A allele of rs6265 
might increase the risk of AD in Caucasian females and female 
LOAD patients. In addition, no evidence indicated an association 
between rs2030324 with AD. Since potential biases and 
confounders could not be ruled out completely in this study, 
further studies are needed to confirm these results. 
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